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A chiral model is introduced that is based on the parity doublet formulation of chiral symmetry 
including hyperonic degrees of freedom. The phase structure of the model is determined. Depending 
on the masses of the chiral partners the transition to the chirally restored phase shows a first-order 
line with critical endpoints as function of chemical potential and temperature in additional to the 
standard liquid-gas phase transition of self-bound nuclear matter. 

We extend the parity doublet model to describe the deconfinement phase transition which is in 
quantitative agreement with lattice data at /j,b = 0. The phase diagram of the model is presented 
which shows a decoupling of chiral symmetry restoration and deconfinement. Loosening the con- 
straint of strangeness conservation we also investigate the phase diagram at net strangeness density. 
We calculate the strangeness per baryon fraction and the baryon strangeness correlation factor, two 
quantities that are sensitive on deconfinement and that can be used to interpret lattice calculations. 

PACS numbers: 21.65.Mn,12.38.Aw,12.39.Fe,25.75.Nq 



I. INTRODUCTION 



The study of dense and hot hadronic matter is a central 
topic of nuclear physics. It is directly linked to the search 
for the phase transition to chirally restored and decon- 
fined matter in ultra-relativistic heavy-ion collisions as 
well as to the study of extremely dense but rather cold 
matter inside compact stars. In spite of several decades of 
experimental and theoretical research the phase structure 
of strongly interacting matter remains uncertain with the 
exception of the regime around cold saturated nuclear 
matter and, to some extent the transition behavior at 
vanishing chemical potential, where lattice gauge calcu- 
lation indicate a cross-over transition to chirally restored 
and deconfined matter, at a temperature currently deter- 
mined to be around 150 to 160 MeV 0,12]. 
At finite chemical potential the phase structure of QCD 
is even less clear. While early extensions of lattice studies 
to finite fiB proposed the existence of a critical endpoint 
at rather small chemical potential [1, Hj], other lattice 
investigations cannot confirm evidence for its existence 



A central point of these investigations is the under- 
standing of the phase transition in the hadronic and 
quark-hadron matter. Recent lattice calculations and 
their analysis in terms of a hadron resonance gas hint to 
the importance of hadronic degrees of matter in driving 
the phase transition to a quark-gluon plasma @, 0, H[. 
Furthermore, the low temperature of the chiral transi- 
tion [2, 12] ! the good agreement with chiral perturbation 
theory below T c [Hi and the apparent sensitivity on the 
hadron properties Q (caused by lattice discretization ef- 
fects) supports the idea that the chiral transition could 
be explained with hadronic interactions. Therefore, also 
a study of purely hadronic models and their properties, 
especially the restoration of chiral symmetry is impor- 
tant. One main benchmark for any useful comprehensive 
model of that kind is a reasonable description of satu- 
rated nuclear matter. In order to have a realistic de- 
scription of highly excited matter strange hadrons have 
to be included in the model description. In a simple linear 
sigma-model it is not possible to have stable bound nu- 
clear matter. Therefore a number of extended approaches 
adding vector and dilaton fields were discussed 
including extensions to flavor SU(3) [l5l - il7| . 
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II. THE PARITY DOUBLET MODEL 

An elegant and alternative description of a transition 
to chirally restored matter is the parity doublet model. In 
this approach an explicit mass term for baryons is pos- 
sible, where the signature for chiral symmetry restora- 
tion is the degeneracy of the baryons and their respec- 
tive parity partners. There are several SU(2) studies of 
nuclear matter adopting this approach showing that it 
is possible to generate saturated matter in the parity 
doublet approach [l8l423j . A SU(3) parity-doublet de- 
scription of hadronic matter was still missing. In (24| 
hyperonic decays in vacuum have been studied in such 
an approach. In the following we outline the basic SU(3) 
parity model. With this ansatz we study nuclear matter 
saturation in order to fulfill one of the benchmarks for a 
useful model as mentioned above. Subsequently we cal- 
culate the phase diagram of isospin-symmetric matter by 
varying the baryonic chemical potential and temperature 
of the system. 

In the parity doublet model positive and negative par- 
ity states of the baryons are grouped in doublets. The 
two components of the fields defining the parity part- 
ners, Lp + and if- transform in opposite way regarding 
chiral transformations: 



ip'_ R = L<p- R (f-L = R^f-L 



(1) 



where L and R are rotations in the left- and right handed 
subspaces. This allows for a chirally invariant mass term 
in the Lagrangian of the general form: 

m (if^j 5 Lp + - ip + j S (P-) = 

m ((f-L(p+R - &-R(p + L - (f+LV-R + (p+RLp-L), (2) 

where mo represents a mass parameter. The general 
SU(3) extension of the approach using the non-linear rep- 
resentation of the fields is quite straightforward as shown 
in [24|]. As outlined in [l5| one constructs SU(3)-invariant 
terms in the Lagrangian including the meson-baryon and 
meson-meson self-interaction terms assuming a nonlin- 
ear realization of chiral symmetry. The part of the La- 
grangian coupling the baryon and the mesonic fields rel- 
evant in a mean-field approximation reads 



Tr(S^H) + m Tr((H 75 T 2 H) +£>WTr(S{S,S}) 
Fi 1 )Tr(S[S,H]) + 5WTr(E)Tr(HH) 
£f>Tr(Hr 3 {£, H}) + F^Tr(ET 3 [E, H]) 
Si 2 )Tr(£)Tr(Sr 3 S) + D v Tr(Z 7 , {V»,Z}) 
F v Tt(E 7ii [V^ S]) + ^Tr(^)Tr(S 7M S) . 



(3) 



Here 5 is the baryon octet whereby each field is a dou- 
blet consisting of the baryon and its negative parity part- 
ner. E and are the multiplets of the scalar and vector 
mesons. The Pauli matrices tv act on the doublets. In 



TABLE I: Model parameters for different values of the mass 
of the nucleonic parity partner 



ko 

(370.99MeV) 2 



fci &2 e 

4.166 -12.902 (75.98MeV) 4 



M N * [MeV] 

gl 
gl 

Ota 

gNu) 



1400 

-7.27 

-0.765 

2.898 

5.375 



1440 

-7.467 

-0.792 

2.812 

5.375 



1490 

-7.714 

-0.823 

2.717 

5.5 



1535 

-7.933 

-0.850 

2.642 

5.563 



general the various sets £)W, F^ l \ £W correspond to the 
D-type and F-type SU(3) invariant baryon-meson cou- 
plings. Note that the parity doublet models allow for two 
different scalar coupling terms i = 1,2. In order not to 
be overwhelmed by coupling constants we will restrict the 
set of non-zero couplings in the actual calculations. As 
the term proportional to mo mixes the upper and lower 
components of the parity doublets, one diagonalizes the 
matrix by introducing new fields B with a diagonal mass 
matrix. Taking along only the diagonal meson contribu- 
tions, the scalar and vector condensates in the mean field 
approximation, the resulting Lagrangian Lb then reads 

i i 

+ ^{Ba^g^+gpiP^ + gti^Bi) . (4) 

i 

The effective masses of the baryons (assuming isospin 
symmetric matter) read 



(g { ai<r + 9$C) 2 + (™o + n s m s y 



0.0 0.1 



(5) 
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FIG. 1: Negative binding energy per particle in MeV as func- 
tion of density of the system for the case Mn* ~ 1535 MeV. 
A reasonable nuclear matter ground state can be achieved. 
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FIG. 2: Phase transition lines of the chiral and liquid gas 
transitions for a mass mjv* = 1535MeV. Solid lines mark 
first-order transitions whereas dashed lines indicate a cross- 
over. The circles mark the critical end points. The liquid- 
gas and chiral symmetry restoration cross-over lines merge at 
higher temperature. 



where the various coupling constants g[^ are given as 
combinations of the original parameters D^\F^\S^ 
in equation [3] and further adding a SU(3) breaking mass 
term that generates an explicit mass corresponding to 
the strangeness n s of the baryon. 

The scalar meson interaction driving the spontaneous 
breaking of chiral symmetry can be written in terms of 
SU(3) invariants h = Tr(£) , h = Tr(£ 2 ) , I 3 = 
det{E) , h = Tr(£ 4 ): 

V = V a + h I 2 - hl 2 2 - k 2 h + fc 3 ln(7 3 ) (6) 

where Vq is fixed by demanding a vanishing potential in 
the vacuum. The explicit symmetry breaking term that 
generates the correct pion and kaon masses with their 
corresponding decay constants can be written as 

L SB = mlU<r + (VZmlfk - -±=m 2 (, 

(7) 

The set of scalar coupling constants are fitted in or- 
der to reproduce the vacuum masses of the nucleon, and 
the Lambda, Sigma, and Xi hyperons, whereas the vec- 
tor couplings are chosen to reproduce reasonable values 
for nuclear ground state properties. The resulting bind- 
ing energy per particle as function of density is shown in 
Fig. 1. For all listed parameters the ground state en- 
ergy per baryon is between -15 MeV and -16 MeV, the 
ground state density has a value of p = 0.15/m^ 3 and 
the compressibility lies between 300 MeV and 310 MeV. 
The latter value is somewhat large. Here, a more detailed 
and extensive parameter study might likely lead to more 
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FIG. 3: Masses of baryons as function of chemical potential. 
The degeneracy of the various parity doublets can be observed 
at high fj,B- 



satisfactory values. Note that the value of the mass pa- 
rameter is set to mo = 810 MeV for all parametrizations. 
Such a choice corresponds to a rather large bare mass of 
the baryons. In principle such a large mo could also be 
generated dynamically through a coupling to the dilaton 
field (see e.g. 25]). Such a coupling can be introduced 
in our model in a straight forward way. In the present 
investigation however we intended to have as few free pa- 
rameters as possible, allthough such an extension will be 
subject of future investigations. 

One candidate for the parity partner of the nucleon is 
the N(1535) resonance. However, this assignment is un- 
clear, the state might also be a broad structure, so essen- 
tially the mass of the particle (assuming its existence) is 
not determined. Resulting parameters for several values 
are shown in Table 1. A SU(3) description, in addition to 
enhance the number of degrees of freedom, also necessar- 
ily increases the number of parameters. In order not to be 
overwhelmed by too many new parameters, for simplicity 
we assume that the splitting of the various baryon species 
and their respective parity partners is of the same value 
for all baryons, which is achieved by setting g^ = g' 2 ' 

and g^ = 0. This should be sufficient for a first inves- 
tigation of the model approach. This assumption agrees 
quite well with the even less certain assignments of the 
parity partners of the hyperons. Obvious candidates are 
the A(1670) and £(1750), whose masses roughly follow 
the equal splitting approximation, assuming the nucle- 
onic parity partner to be the N(1535). In the case of the 
S* the data are unclear. 

In another simplification the hyperonic vector inter- 
actions were tuned to generated reasonable optical po- 
tentials of the hyperons in ground state nuclear matter, 
with U A {p ) -28 MeV and U s (p Q ) - 18 MeV . The value 
for the strange quark mass was fixed at m s = 150 MeV. 
The numbers used are summarized in Table 1. A more 
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FIG. 4: (Color online) Normalized value of the chiral conden- 
sate as a function of temperature, at fin = 0. The solid line 
is the model result for To = 220 MeV and the dashed line for 
To — 270 MeV. The grey dashed line depicts results for the 
pion contribution to the chiral condensate from chiral pertur- 
bation theory [2(|. The symbols denote lattice data, where 
the open symbols represent older data with the asqtad and p4 
action and the colored symbols are more recent results (see 
text). 



FIG. 5: (Color online) Value of the Polyakov loop as a func- 
tion of temperature at /xs =0. The solid line is the model 
result for T = 220 MeV and the dashed line for T = 270 
MeV. The symbols denote lattice data, where the open sym- 
bols represent older data with the asqtad and p4 action and 
the colored symbols are more recent results. 



III. INCLUDING MESONS AND QUARKS 



exhaustive study of various parameter setups will be per- 
formed in future work. 

The equations of motion following from Eqs. (4, 6, 
7) are then solved self-consistently in mean field approx- 
imation by minimizing the grand canonical potential as 
function of baryonic chemical potential and temperature. 

The resulting phase diagram for the transition from 
chirally broken to chirally restored phase is shown in Fig. 
2. One can observe two first-order transition lines at high 
densities. The first one is the liquid-gas phase transition, 
which indicates that the model exhibits a bound nuclear 
ground state. The second one is the chiral transition, 
also signaling the onset of the population of the baryonic 
parity partners. The lines stop at second-order critical 
end points, at (T c ,/i c ) = (21 MeV, 1560 MeV) and (17.5 
MeV, 905 MeV) . The critical point of the chiral transition 
is very low in temperature. At values of Mm* below 1460 
MeV the first-order transition becomes a cross-over for 
all values of T and fi. Both cross-over lines join at a 
temperature T « 120 MeV. 

Fig. 3 shows the effective masses of the baryons with 
their parity partners as function of baryochemical poten- 
tial. One can observe the effect of the two phase tran- 
sitions, leading to essentially degenerate opposite parity 
states. 

Here the calculations were done for isospin symmetric 
matter with strangeness zero, which would be the logi- 
cal first-order assumption for matter created in heavy- 
ion collisions. The study of star matter in beta equilib- 
rium including leptons for ensuring charge neutrality is 
in progress. Results will be presented in a forthcoming 
publication. 



Recent results from different lattice collaborations in- 
dicate that the chiral phase transition at [i q = occurs at 
a rather low temperature. Although there are still con- 
siderable systematic lattice effects there seems to be a 
consensus that the crossover temperature is between 150 
and 160 MeV. This can be seen in figure (|4]) where the 
symbols denote the different lattice actions. In particu- 
lar intriguing is the fact that the value of the chiral con- 
densate already drops to about 80 to 90% of its vaccum 
value at temperatures less than 150 MeV. In the parity 
doublet model such an early decrease of the chiral con- 
densate can be hardly accommodated by only coupling 
the baryons to the chiral fields as their mass is rather 
large and they are not thermodynamically activated at 
such low temperatures. On the other hand calculations 
in chiral perturbation theory have shown that the pionic 
contribution to the chiral condensate is considerable at 
such low temperatures because of the small pion mass. 
This is also shown in figure (U]), as the grey dashed line 
depicts the chiral perturbation theory results for the pion 
self interaction taken from (26j. One can see, that at low 
temperature the behavior of the chiral condensate seems 
dominated by the pseudoscalar contributions while only 
at larger temperature the baryon interactions become im- 
portant. Consequently we would like to include effects of 
the pseudoscalars in our parity doublet model to accom- 
modate the low temperature behavior of the chiral con- 
densate at small net baryon densities. To this end we take 
into account the coupling of the scalar field, which origi- 
nates from the explicit symmetry breaking term, Eq. (7) 
(see Ref. [27|). Including the pseudoscalar mesons this 
term generates a mass for the pions as = rn^ a/ao. 
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FIG. 6: (Color online) The interaction measure (denned as 
(e— 3p)/T 4 ) as a function of temperature at /is = 0. The solid 
line is the model result for To = 220 MeV and the dashed line 
for To = 270 MeV. The symbols denote lattice data, where 
the open symbols represent older data with the asqtad and 
p4 action and the colored symbols are more recent results. 
The grey line represents the interaction measure for a hadron 
resonance gas equation of state which includes all hadronic 
resonances up to 2.2 GeV. 



This leads to an increase of the pressure of the pion gas 
with decreasing scalar condensate, thus driving the phase 
transition to lower temperatures. 

It is well known that at some temperature QCD ex- 
hibits a transition to a deconfined phase at which the 
quarks become the dominant degrees of freedom. When 
this deconfinement will appear and what the order pa- 
rameter for this transition might be is still under heavy 
debate [H, [2{| . Assuredly one can only say that it oc- 
curs in a temperature region of Td ec ~ 160 — 400 MeV. 
Nevertheless at some point the hadronic parity doublet 
model will not be the appropriate effective description of 
QCD and one needs to introduce a deconfinement mech- 
anism in the model. In this work we will apply a mech- 
anism that has been introduced in [3(| to add a decon- 
finement transition in a chiral hadronic model. This is 
done by adding an effective quark and gluon contribu- 
tion as done in the PNJL approach py], [32j ■ This model 
uses the Polyakov loop $ as the order parameter for de- 
confinement. $ is defined via $ = |Tr[exp(z J drA^)], 
where A4 = iAq is the temporal component of the SU(3) 
gauge field, distinguishin g <i> and its conjugate $* at fi- 
nite baryon densities [33-35]. In recent years the PNJL 
model has been widely used and extended to include non- 
local interactions as well as an imaginary chemical poten- 
tial (see also [36|-[65j]). 

The effective masses of the quarks are generated by 
the scalar mesons except for a small explicit mass term 
(8m q = 5 MeV and 8m s = 150 MeV for the strange 



FIG. 7: (Color online) Number densities of the different par- 
ticle species as a function of temperature at /is = 0. The 
solid lines show results for To = 220 MeV and the dashed 
line for To = 270 MeV. The mesons are shown in green, the 
densities for baryons plus anti-baryons in red and the quarks 
plus anti-quarks in black. 

quark) and mo (explain): 

m * q = 9q<j<J + 5m q + m 0q , 

m l = 9s(( + 8m s + m 0q , (8) 

with values of g qa = g s ^ = 4.0. As in the case of the 
baryons we also introduced a mass parameter mo q = 200 
MeV for the quarks. Again this additional mass term 
can be due to a coupling of the quarks to the dilaton 
field (gluon condensate). Fot this mass term the quarks 
do not appear in the nuclear ground state which would 
be an unphysical result. This allows to set the vector 
type repulsive interaction strength of the quarks to zero. 
A non-zero vector interaction strength would lead to a 
massive deviation of the quark number susceptibilities to 
lattice data as has been indicated in different mean field 
studies [gi-ii]. 

A coupling of the quarks to the Polyakov loop is intro- 
duced in the thermal energy of the quarks. Their thermal 
contribution to the grand canonical potential il, can be 
written as: 

^ = -Tg^/^ln(l + $exp^i) (9) 
and 

^=- T Ej^J ^ In (l + exp 

The sums run over all quark flavors, where 7$ is the 
corresponding degeneracy factor, E* — m* 2 + p 2 the 
energy and /z* the chemical potential of the quark. 
All thermodynamical quantities, energy density e, en- 
tropy density s as well as the densities of the different 
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particle species pi, are derived from the grand canonical 
potential. It includes the effective potential £/($, $*, T), 
which controls the dynamics of the Polyakov-loop. In our 
approach we adopt the ansatz proposed in (32j : 

2 v ; 

+ 6(T)/n[l-6$$*+4($ 3 $* 3 )-3($$*) 2 ] (11) 



feT 3 T. 



with o(T) = a Q T 4 + ai T T 3 + a 2 T$T 2 , b(T) = u 3 j 
The parameters <Zo,ai,a2 and 63 are fixed, as in [32j . 
by demanding a first order phase transition in the 
pure gauge sector at T = 270 MeV, and that the 
Stefan-Boltzmann limit of a gas of gluons is reached for 
T — > 00. In general of course the presence of quarks may 
have a significant influence on the Polyakov potential 
[42} so one should not regard the parameters to be 
absolutely fixed. 

In the following we introduce excluded volumes 
for the hadrons in the system. As a onsequence the 
hadronic contributions from the equation of state at high 
temperatures and densities will be suppressed. Including 
effects of finite-volume particles in a thermodynamic 
model for hadronic matter, was proposed long ago 
[69Mrsl ] . In recent publications [3(3, HI[ we adopted this 
ansatz to successfully describe a smooth transition from 
a hadronic to a quark dominated system (see also 79]). 

In particular we introduce the quantity Vi which is 
the volume excluded of a particle of species i where we 
only distinguish between hadronic baryons, mesons and 
quarks. Consequently Vi can assume three values: 

VQuark = 
v Baryon v 
VMeson = v/o, 
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FIG. 8: (Color online) Densities of the different particle 
species at T = and as a function of fi q — [ib/3. The black 
solid line depicts the density of protons plus neutrons, the 
red dashed line for up plus down quarks and the green dash 
dotted line for the chiral partners of the nucleons. 



= f e ' 1 Pi = f Pi s = ^f 



(14) 



As a consequence, the chemical potentials of the 
hadrons are decreased by the quarks, but not vice 
versa. In other words as the quarks start appearing 
they effectively suppress the hadrons by changing their 
chemical potential, while the quarks are only affected 
through the volume correction factor /. 



IV. THERMODYNAMIC PROPERTIES 



Where a is a number larger than one. In our cal- 
culations we choose a value of a = 8, which assumes 
that the radius r of a meson is half of the radius of a 
baryon. Note that at this point we neglect any possible 
density-dependent and Lorentz contraction effects on the 
excluded volumes as introduced in [jf| [73] ■ 

The modified chemical potential Jli , which is connected 
to the real chemical potential pi of the z-th particle 
species, is obtained by the following relation: 

Pi = Pi ~ v iP (12) 

where P is the sum over all partial pressures. To be 
thermodynamically consistent, all densities (e^, pi and 
Si) have to be multiplied by a volume correction factor /, 
which is the ratio of the total volume V and the reduced 
volume V', not being occupied: 

/ = v = (i+E^) _1 ( 13 ) 



A surprising result from recent lattice studies of the 
2+1 flavor QCD equation of state at finite temperature 
is the apparent decoupling of the chiral phase transition 
and the increase of the Polyakov loop which was long 
thought to be a good order parameter for deconfinement. 
In particular one observes that the steepest change in the 
chiral condensate occurs at a low temperature of roughly 
150— 160 MeV, depending on the choice of the lattice ac- 
tion as well as the scale which is used to translate lattice 
quantities into the physical temperature. On the other 
hand a considerable increase in the value of the Polyakov 
Loop is only observed above T ss 210 MeV, where this 
result seems quite independent of the lattice action that 
is applied. 

In the following we will compare results on the or- 
der parameters and thermodynamic quantities calculated 
within our model to recent lattice results from the differ- 
ent collaborations [3, H i, E3, IMQ Figure U displays 
the results for the temperature dependence of the expec- 
tation value of the a field, normalized to its ground state 



7 



T =270 MeV 



400 



300 



CD 



200 



100 




400 400 



100 200 300 400 500 

u [MeV] 



T =220 MeV 

100 200 300 400 50 



300 300 



200 200 



100 100 



41 




00 



300 



200 



100 



CD 



100 200 300 400 500 

]i [MeV] 



FIG. 9: (Color online) Phase diagrams for our model with two different values of To. Within the orange region the normalized 
value of the chiral condensate lies between 0.2 < cr/co < 0.8. In the green area the value of the Polyakov loop is in between 
0.2 < <f> < 0.7. The black lines indicate first order phase transitions, where the points indicate the critical endpoints. 



value, from the SU(3) parity model including quarks. 
The dashed line depicts the results when the parameter 
of the Polyakov potential is unchanged (To — 270 MeV) 
and the solid line when T is changed to 220 MeV. Note 
that such readjustments of To are commonly used in a 
number of PNJL studies and To can even depend on \ib 

Our results are compared to calculations from chiral 
perturbation theory (grey short dashed line) [26| and 
recent lattice results (colored symbols depict the more 
recent results while the open symbols refer to previously 
used lattice actions). 

At temperatures below 160 MeV the decrease of the 
a is dominantly caused by the pseudoscalar-scalar 
coupling. Compared to the chiral perturbation theory 
results, which depict the pion contribution to the chiral 
condensate, our model still shows a slower decrease of 
the chiral condensate with temperature, while lattice 
data are rather well described by the chiral perturbation 
theory up to temperatures of 150 MeV. Above that 
temperature the baryon scalar interaction (and at even 
higher temperature the quark scalar interaction) starts 
to contribute to the value of the cr-field. In this regime 
our model gives a good description of the lattice data. 
The chiral critical temperature T c x only weakly depends 
on the value of To because chiral symmetry restoration 
is mainly driven by the hadronic interactions with the 
fields. We obtain the values of T c x from the maximum in 
da/dT as T* = 172 and 165 MeV respectively. 



While the temperature dependence of the lattice 
results on the chiral condensate strongly depends on 
the action and lattice spacing that is applied, this 
dependence is not so strong for the Polyakov loop as can 
be seen in figure [5] Here again we compare our model 
results with lattice data. Because To directly influences 
the dynamics of the Polyakov loop we see a strong 
dependence of T C PL on the change in To. We obtain 
T c pi = 210 and 173 MeV. the general shape of the curve 
of 4>(T) is not changed considerably with To and we 
observe that we cannot accommodate the almost linear 
increase of the Polyakov loop at lower temperatures, as 
seen on the lattice, by simply adjusting Tq. 



Recent lattice studies start to agree on the tempera- 
ture dependence of the order parameters. However this 
is not the case when the interaction measures, closely 
related to the thermodynamic properties of the matter, 
are compared. In figure [6] we compare our model results 
on the interaction measure, for the two values of Tq, with 
different lattice data sets. First we have to note that 
the results from the HotQCD collaboration (red squares 
and green triangles) differ considerably from those 
of the Wuppertal-Budapest group (black diamonds). 
Furthermore we see that the interaction measure for 
our calculation with To = 220 MeV gives a drastically 
better agreement with the HotQCD lattice results than 
for 270 MeV. in any case our model underestimates 
the contribution to the interaction measure at low 
temperatures. However, this could be understood by 
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FIG. 10: (Color online) Value of the Polyakov loop as a func- 
tion of temperature at \ib = 0. This time the solid line is 
the model result for To = 210 MeV and a modifies parameter 
a\ — —8.47. The dashed line stands for the previous result 
with T = 270 MeV and oi = -2.47. The symbols denote lat- 
tice data, where the open symbols represent older data with 
the asqtad and p4 action and the colored symbols are more 
recent results. 



missing contributions from hadronic resonances as is 
seen from the grey line which depicts (e — 3p)/T for 
a hadron resonance gas which has already been shown 
to give good results for the interaction measure at 
temperatures below T*. 

The transition from a hadron to quark dominated system 
is depicted in figure [7] with the different particle number 
densities at hb =0. While the mesons dominate the low 
temperature region (green lines) one can clearly see that 
the hadrons are slowly removed from the system as the 
quark density increases (black lines) and consequently 
the hadrons play less of a role for the thermodynamic 
quantities. 



A. Finite \ib 

An advantage of our effective model is that, unlike 
lattice studies, we can easily extend our studies to fi- 
nite baryon densities and explore the phase behavior at 
[iB > 0. Figure [5] shows again the number densities of 
different hadrons and quarks, this time at T = 0, and 
as a function of \xb- at the liquid-gas phase transition 
the nucleon density exhibits a jump as expected. At 
even higher chemical potentials we can observe a sec- 
ond jump in the densities, at this point the chiral part- 
ner of the nucleon is activated, as well as the quark de- 
grees of freedom. Both steps in the densities correspond 
to jumps in the chiral condensate as can be seen more 
clearly in figure [SJ In this figure we depict the phase 
diagrams obtained from our model for the two values 
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FIG. 11: (Color online) The interaction measure (defined as 
(e — 3p)/T 4 ) as a function of temperature at fis = 0. This 
time the solid line is the model result for To = 210 MeV and 
a modifies parameter ai = —8.47. The dashed line stands for 
the previous result with T = 270 MeV and ai = -2.47. The 
symbols denote lattice data, where the open symbols repre- 
sent older data with the asqtad and p4 action and the colored 
symbols are more recent results. The grey line represents the 
interaction measure for a hadron resonance gas equation of 
state which includes all hadronic resonances up to 2.2 GeV. 



of To. The black lines with endpoints depict the region 
where the change of the chiral condensate is of first or- 
der. Instead of drawing an ambiguous crossover 'line' we 
show the regions in which the value of the chiral con- 
densate changes from 0.2 < cr/co < 0.8 as the orange 
area. This region we will refer to as the chiral crossover 
region. The green area is defined as the region where the 
Polyakov loop is between 0.2 < cf) < 0.8, which we will 
refer to as the deconfinement crossover. Please note that 
within our model the Polyakov loop generally changes in 
a smooth crossover and only exhibits a first order phase 
transition at the chiral phase transition, where it jumps 
from cj> — to 0.1. The position of the chiral critical 
endpoint is T cep = 58 MeV and f200 MeV for 

T = 270 MeV and for T = 220 MeV only the chemical 
potential changes to /j,g P « 1150 MeV. Such a low tem- 
perature endpoint as well as the fact that the crossover 
region does not become considerably wider with increas- 
ing chemical potential both agree well with recent lattice 
findings Q. Generally the deconfinement crossover oc- 
curs at a higher temperature (higher chemical potential) 
as the chiral phase transition. This gives rise to an in- 
teresting phase of chirally symmetric confined hadronic 
matter. 



B. Varying Polyakov Loop parameters 

As had been mention in the previous section, our model 
does not yield a good description of the temperature be- 
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FIG. 12: (Color online) Contour plots of the normalized chiral condensates as a function of the chemical potentials [ib and /is 
for fixed temperature. The red lines correspond to different values of a fixed strangeness to baryon fraction /„. 



havior of the Polyakov loop together with a good descrip- 
tion of the interaction measure. Within the PNJL model 
it has also been pointed out that a concurrent description 
of both order parameters and the interaction measure is 
usually not achieved when simply To is varied (see e.g. 

In our model in particular the linear increase of the 
Polyakov loop is ill described. Usually the parameters 
of the Polyakov potential are fitted to pure gauge lattice 
result. However it has been shown that the presence of 
quarks may have an influence on these parameters (42j 
and in general there is no reason that only To should be 
affected by such a quark coupling. To investigate and 
illustrate what effects a different parameter set of the 
Polyakov potential has on the interaction measure we 
will adjust one parameter (ai) to ai = —8.47 (instead 
of ai = -2.47). Together with changing T to T = 210 
MeV this gives an improved description of the lattice 
data for the Polyakov loop behavior as can bee seen in 
figure [TUl 

Figure [11] shows the resulting interaction measure for our 
adjusted parameter. One can see that, due to the slower 
rise of the Polyakov loop, the peak in the interaction mea- 
sure is considerably lowered. Taking into account miss- 
ing contributions from resonances one can even conclude 
that this parametrization compares more favorably with 
the Wuppertal-Budapest results than with the HotQCD 
lattice data. 

Consequently, to understand the interplay between the 
order parameters and the interaction measure from lat- 
tice calculations, one has to resolve the still existing dis- 
crepancies. Only then conclusions regarding the role of 
the degrees of freedom in the chiral and deconfmement 
phase transition can be drawn. 



V. THE STRANGE PHASE DIAGRAM 

Until now all calculations where restricted to the limit 
of vanishing net strange density. Because the strong in- 
teraction conserves the net strange particle number, the 
equation of state used for the description of heavy ion 
collisions is usually considered to be net strange free. 
However, there are several issues that make the study 
of the strange EoS interesting. Some of these issues are: 

1. As has been shown in [87j the net strangeness dis- 
tribution in coordinate and momentum space of a 
heavy ion collision can fluctuate, although the to- 
tal net strangeness is zero. To dynamically treat 
such a system, and calculate observables that arise 
from such a strangeness fluctuation, the equation 
of state for p s ^ needs to be evaluated. 

2. Compact stars are very dense and long lived ob- 
jects. Due to a /^-equilibrium inside the star, net- 
strange conservation can be violated by the weak 
interaction. 

3. Lattice QCD results at finite /is are often evaluated 
through a Taylor expansion in Hb at Hb = Ms = 
0. A vanishing strange number chemical potential 
usually induces a non- vanishing net strangeness, 
which means that the equation of state of net- 
strange matter is calculated. 

First investigations on the strange equation of state 
were done in |8a |. where one usually considered a first 
order transition from a hadron to a quark phase. In 
our model we are able to discuss the strange EoS in the 
context of a smooth transition from a confined hadron 
phase to a deconfined quark phase. 

Figure [T2] presents our results on the order parameter 
of the chiral phase transition as a function of hb and 
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FIG. 13: (Color online) Baryon number density divided by 
T s as a function of temperature for hb /T = 1. Shown is the 
exact solution from the model (black solid line) and Taylor 
expansions of the density, taking into account only the second 
(green dashed line) and forth (red solid line) coefficient. 



/is at fixed temperature. The red lines indicate paths 
of constant values for f s = p s /pB, the strangeness per 
baryon fraction, note that f s = corresponds to our 
results in section [TV] (with T = 220 MeV). At the 
temperature T = 56 MeV, the critical endpoint of the 
chiral phase transition was located at « 1150 MeV. 
We can observe that for increasing f s , the change in the 
order parameter becomes steeper and the value of Tcep 
increases slightly to T CE p = 68 MeV for f a = 0.5. At 
the larger temperature we also observe a slight change 
in the phase structure. Here, for increasing f s , the 
crossover moves closer to the ^b,s = 0-line. 

For a gas of deconfined quarks there is a strong corre- 
lation between the baryon number and strangeness. In a 
hadronic medium such a correlation is usually not triv- 
ial as strangeness can be found in mesons and baryons. 
These considerations led to the idea that the so called 
strangeness-baryon correlation factor cbs is sensitive to 
the deconfinement and/or chiral phase transition [89| . 
On the other hand the strangeness to baryon ratio f s 
should also be sensitive on any phase transition at finite 
baryon densities. 

On the lattice such quantities are usually calculated as 
functions of the expansion coefficients, it is defined as 



cbs 



, (N B N S ) - (Nb) (N s ) 
(Nl) - (N S ) 2 



(15) 



The question is how many coefficients are needed to 
evaluate the baryon and/or strange densities at finite 
Pb/T to a given accuracy [§oj]. As an example figure HH 
shows the baryon density as a function of temperature 
at fixed hb/T = 1 and ps = for our exact model 
calculation. Alternatively we can also numerically 



FIG. 14: (Color online) Shown are the strangeness to baryon 
correlation coefficient cbs (red short dashed line) compared 
to the quark-gluon fraction A = eQ uarfes+G;ilona /e To t (black 
solid line) as a function of temperature for (ib/T — 3 and 
[is = 0. The plot also shows the strangeness per baryon 
fraction f s (green dashed line) and the quark-gluon fraction 
A (grey solid line) as a function of temperature for us /T — 1 
and hb = 0. 



extract the expansion coefficients for our model and 
expand the density in powers of [ijT . One can see that 
already the result, taking into account only the first 
non vanishing coefficient C2, gives a quite reasonable 
description of the exact result. Taking into account the 
2nd and 4th order coefficient already allows to describe 
the exact result to high accuracy, except at the point of 
the crossover transition. This means that, in order to 
calculate f s at finite fJ-s/T, it is sufficient to extract the 
coefficients up to 4th order (c^'f) from the lattice. 

The information that can be extracted from these 
quantities is exemplified in figure 1141 Here we show the 
exact solution for cbs a $ a function of temperature for 
Pb/T = 3 and ps = 0. One can observe a distinct 
peak at T w 150 MeV => fi B = 450 MeV. Comparing 
with figures [9] and [12] one can identify this peak with 
the crossover transition of the chiral condensate. Such 
a behavior of cbs has been predicted and also has been 
shown to exist in lattice data 



9l|. At higher temper- 



atures the strangeness to baryon correlation approaches 
unity which resembles closely the behavior of the quark 
and gluon fraction A = eQ uar ks+Giuons / &Tot of the sys- 
tem. In comparison figure 1141 also shows the tempera- 
ture dependence of f s at ps/T = 1 and [Lb — 0. This 
quantity is even more sensitive in the quark-gluon frac- 
tion as is cbs, while it seems to be not very sensitive to 
the chiral phase transition. The peak in f s can rather 
be understood as a consequence of our excluded volume 
treatment, where mesons have a smaller excluded volume 
than baryons. Hence mesons, that can carry strangeness, 
are less suppressed than baryons. 
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VI. DISCUSSION 

We presented results on the phase structure of a SU(3) 
parity-doublet description of hot and dense hadronic 
matter. With appropriate parameters we could generate 
a quantitatively acceptable nuclear ground state. The 
phase diagram in temperature and baryochemical poten- 
tial exhibits a liquid-gas first-order phase transition as 
well as a chiral phase transition that is connected to the 
population of the parity partners and the onset of their 
degeneracy with the normal baryon states. Depending 
on the mass gap between the baryons and their parity 
partners this transition is first-order at high densities 
and low temperatures and a crossover otherwise, or a 
smooth crossover for all values of T, /i for smaller mass 
gaps. In order not to be overwhelmed by too many new 
parameters some simplifications of the parameter choice 
have been made, assuming an equal mass gap between all 
positive and negative parity baryons. These restrictions 
should be relaxed in further studies to explore the model 
in more detail. 

In the second part of this paper we extended the SU(3) 
parity doublet model to incorporate a deconfinement 
phase transition. When comparing our results to lattice 
data at fiB — we find that the low temperature 
behavior of the chiral condensate is dominated by 
hadronic interactions. In such a scenario a decoupling of 
the Polyakov loop and the chiral condensate, as is seen 
in recent lattice studies, can be easily understood. A 
feature which is common to PNJL-type models is that 
a simultaneous description of the interaction measure 
and the Polyakov loop cannot be achieved simply by 
adjusting the parameter Tq. If we loosen also the 
constraint on other parameters of the Polyakov loop 
potential we obtain an improved description of the 
Polyakov loop dynamics. It would be interesting to 
investigate if, e.g. in PNJL models, the slower increase 
also shows to have drastic effects on the interaction 
measure as has the parameter change presented in this 
work. As lattice results still differ strongly in their 
results on the interaction measure it is not possible to 
say if such a reparametrization improves or weakens 
the model. Consequently it is of utmost importance 
to understand and settle the differences in the lattice 
results to be able to understand the interplay between 
the order parameters and the thermodynamics, i.e. the 
active degrees of freedom. 

At finite baryon densities our model describes the 
deconfinement transition as a continuous crossover 
for all values of hb- Only the chiral order parameter 



exhibits two discontinuities. One is related to the 
nuclear liquid gas phase transition while the other can 
be identified as the chiral phase transition and appears 
at larger densities. We also observe that the critical 
endpoint of the chiral phase transition has a rather 
small temperature T cep = 56 MeV. As the chiral phase 
transition is driven mainly by hadronic interactions 
and the deconfinement by quarks and the Polyakov 
potential, we see a decoupling of the order parameters, 
which becomes stronger for large chemical potentials. 
We observe several different states of matter that can 
form, starting from a nucleon liquid which changes to 
a phase of chirally symmetric hadrons. Only at higher 
temperature these hadrons disappear and the quarks 
are the dominant degrees of freedom. Whether such a 
chirally symmetric hadronic phase can be the A^ c = 3 
equivalent of the N c = oo quar kyonic phase [92[ is still 
under extensive debate |93l - l95l ]. In any case the high 
density part of the QCD phase diagram could have a 
rather rich phase structure to explore. 

In the last part of this paper we discuss the properties 
of our models phase diagram at finite net-strange den- 
sity. This aspect of QCD matter is not only interesting 
for heavy ion collisions and compact stars, but also for 
a comparison with lattice results extrapolated to finite 
M-B.s- We find that the location of the critical endpoint 
shifts to a slightly higher temperature for a finite net 
strangeness (lattice results). 

We briefly discussed quantities that are sensitive on the 
chiral and/or deconfinement phase transition. In partic- 
ular these are the strangeness baryon correlation factor 
cbs an d the strangeness per baryon fraction f s . Both 
show to be sensitive to the deconfined fraction on the 
system while cbs also shows a distinct peak at the chi- 
ral crossover at finite chemical potential. The advantage 
of extracting / s from the lattice is that it can be evalu- 
ated e.g. in a Taylor expansion, using only the first two 
non-zero expansion coefficients for the strange and light 
quark number susceptibilities. 
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